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Abstract: 

In this article we study effects of parity nonconservation in atomic dysprosium, where one 
has a pair of nearly degenerate levels of opposite parity. We consider the time evolution of 
this two-level system within oscillatory electric and magnetic fields. These are chosen to have 
a periodical structure with the same period, such that a Floquet matrix describes the time 
evolution of the quantum states. We show that, if the states are unstable, the eigenvalues 
of the Floquet matrix may have contributions proportional to the square root of the parity 
violating interaction matrix element H^ while they are almost degenerate in their parity 
even part. This leads to beat frequencies proportional to \/H„ which are expected to be 
larger by several orders of magnitude compared to ordinary P-violating contributions which 
are of order H„. However, for the simple field configurations we considered, it still seems 
to be difficult to observe these P-violating beat effects, since the states decay too fast. On 
the other hand, we found that, within only a few Floquet cycles, very large parity violating 
asymmetries with respect to experimental setups of opposite chirality may be obtained. The 
electric and magnetic fields as well as the time intervals necessary for this are in an experi- 
mentally accessible range. For statistically significant effects beyond one standard deviation 
a number of about 10'' atoms is required. Our ideas may be applied directly to other 2-level 
atomic systems and different field configurations. We hope that these ideas will stimulate 
experimental work in this direction. 
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1 Introduction 



In this article we study theoretically possible ways to obtain large parity (P) violating 
effects in a time-dependent two-level quantum mechanical system. This system is 
described by a Hamiltonian operator of a fairly general form. However its choice has 
been motivated by studies of P- violating effects in atomic dysprosium [|l] . 

Nowadays, using different elements (Cs, Tl, Pb, Bi), atomic parity non- conservation 
(PNC) measurements have already reached high precisions, for the case of Cs of ~ 0.3% 
experimentally, with up to 1% theoretical uncertainty [0. (For reviews of PNC in 
atomic physics cf. |^.) Dysprosium is another interesting system where atomic PNC 
can be measured ([0], [^). There is an enhancement of the PNC mixing of a pair of 
opposite parity levels due to their near degeneracy. The two levels (4/^° 5d6s, J = 10), 
and (4/^5(i^6s, J = 10) (cf. Fig. 1 of [Q]) are split by an amount of the magnitude 
of the hyperfine splittings. In particular, the hyperfine components of ^^^Dy (nuclear 
spin: I = 5/2) with total angular momentum F = 10.5, the closest pair of opposite 
parity levels, are split by only A/h = {Eb — EjCj/h = 3.1 MHz while the separation of 
them to other levels is at least of the order of GHz (cf. Fig. 2 of |l|]). Moreover with 
dysprosium one may reduce the experimental uncertainties by comparing the effect for 
its different stable isotopes. 

Consider now these states in external electric and magnetic fields £ = Se^, 3 = Be^, 
respectively. Here is the unit vector in 3-direction, chosen as the quantization axis, 
and S and B are allowed to vary in time. For suitable values of the magnetic field 
B, the sublevels of the opposite parity states with F = 10.5 and a certain magnetic 
quantum number, e.g. mp = —10.5, are nearly degenerate, i.e. the energy separation 
between them is much smaller than the separation between either of these levels and 
any other sublevel. Then we may exclusively consider this 2-level subspace TZ with 
basis 

\A) = |4/io 5rf 6s, F = 10.5, = -10.5), 

\B) = |4/9 5c/2 6s, F = 10.5, = -10.5). (1.1) 

The dynamics in the subspace TZ is governed by an effective Hamiltonian operator, 
which is represented by a non-hermitean 2x2- matrix 

TT ^(-^iaB-^Ta d£ + iH^ \ 

-^ff \ d£-iH^ A-fiBB-iTBj' ^ ' 

where the first row and column refers to level A, the second to B. The effective 
Schrodinger equation 

d_ 

di 

describes the time evolution of the projection \t) of the total time dependent state 
vector onto the subspace TZ within an approximation where the states \A) and \B) 
decay exponentially. The corresponding decay widths are T a,b- With the standard 
choice of phases for the states (|1 . 1| ) , the P- violating matrix elements 

{A\H^^o\B) = -{B\H^^c\A) = iH^ (1.4) 



^-Jt) = H,ff\t) (1.3) 
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are purely imaginary due to time reversal (T) invariance. The real constants /ia and 
liB are related to the electric dipole moment connecting the levels A and B and the 
magnetic dipole moments respectively: 

d{mF) ■■= -{B,F,mF\D-i\A,F,mF) = -iA,F,mF\D3\B,F,mF) 

= V!^i£l^^l^4£^(5||Z)||A), (1.5) 
Jj{J + 1){2J +1) 



liaijnp) '■= {a, F,mF\fi3\(y, FjUIf) 

eh 

= ~9a,FmF^^, {ae{A,B}). (1.6) 

The ( J, /; mj, mj\F, tuf ) are Clebsch-Gordan coefficients. The reduced matrix element 
of D was measured to be | (i?| |D| |y4)| = 1.5(1) ■10~^e as |Q, where is the Bohr radius 
and — e and rrie are the charge of the electron and its mass, respectively. The Lande 
factors ga,F are given in terms of the factors ga,j @], ga,F = Qa,.] [F{F + 1) + J {J + 
1) - /(/ + 1)]/[2F(F + 1)], which in turn, are gA,j = 1-21, gsj = 1-367 (cf. [| and 
P]). Numerical values for the constants, which we will use in this paper, are given as 
follows (cf. i): 

TA/h = 20 kHz, 
TB/h = IkHz, 
d{-10.5)/h = 3.8 kHz V-^cm, 
A/h = 3.1MHz, 
l2Ai-lO-5)/h = 16.4MHz 
/iB(-10.5)//i = 18.5MHz G-\ 

H^/h = 0.002 kHz. (1.7) 

Note that in |jl|, Ts/h < IkHz was given.0 The value for represents roughly the 
result for the upper limit on H„ reported by |T| (cf. their Fig. 14). We assume this 
value of H„ for illustration. Changing to other values of H„ is straightforward. 



One common procedure to measure PNC in atomic systems is the following: One 
prepares an atomic beam with the atoms being in some definite initial state. This state 
evolves in time as the beam passes through some electromagnetic field configuration. 
Occupation numbers corresponding to the projection of the resulting final state onto 
some subspace are measured. If the whole setup has handedness, i.e. is non-invariant 
under a certain improper rotation, PNC effects manifest themselves as changes of such 
occupation numbers under a reversal of the handedness of the setup. 

The observables in such an experiment are thus occupation numbers of states, which 
should have contributions which change under a P transformation if PNC effects are to 
be found. For instance, such contributions could be, to leading order, linear in the P vi- 
olating matrix element For the above case of dysprosium the time evolution of the 

^ In past experiments, the lifetime of the state showed to be longer than the time of flight of an 
atom within the apparatus . 
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states ( |1.1D was studied in []T| for a rapidly varying electric field, keeping the magnetic 
field constant. The frequency of the electric field was chosen large enough such that 
the changes in the system were non-adiabatic. From their measurements they report 
an upper limit on the P-violating matrix element of \H„\/h = \2.4:±2.9stat ^^-^systl Hz 
Note that the present theoretical prediction is H^/h = 70 ± 30 Hz 0. 



The goal of the present paper is to investigate possible enhancement mechanisms 
for PNC effects in dysprosium. Extensions to other systems are straightforward. For 
a similar experimental configuration as it as been used in fl]] we consider observables, 
whose sensitivities are large enough to show large P-violating effects and might allow 
for measurements of the PNC parameters with a higher precision. The underlying 
ideas have been developed within earlier work on P-violating polarization rotations 



and P-violating energy shifts in light hydrogen-like atomic systems 0, ^ |Tl|]. The 
main new idea now is to study the Floquet operator, i.e. the time evolution operator 
for one cycle of an external potential being periodic in time. It is well known that, 
knowing the eigenvalues and eigenstates of this operator and the time evolution within 
one period, it is in general possible to calculate the evolution for all times [0 . In order 
to simplify our calculations we will consider piecewise constant electric and magnetic 
fields, which vary discontinuously in time, and use the sudden approximation for the 
transitions to compute the time evolution of the atomic states. We will choose the 
fields such that the if„-independent parts of the Floquet-eigenvalues are degenerate 
and investigate the possibility to observe PNC effects which are proportional to the 
square root of H^. The result will be, that, for specific initial states, it is possible to 
observe, after a definite number of Floquet-periods, very large asymmetries, of order 
100%, in the occupation number of one of the states 1^4), \B) w.r.t. the reversal of the 
fields' handedness. 



The organization of our paper is as follows: In Section || we outline the formalism 
for the description of the time evolution of the atomic states using the Floquet oper- 
ator. In Section ^ we calculate the eigenvalues and their P-conserving and -violating 
contributions for a specific electromagnetic field configuration and study the role of 
the non-hermitecity of the Hamiltonian. In Section we present our PNC observables 
and consider their statistical significance. Finally, in Section |^, we present numerical 
results for the above mentioned states of dysprosium and give our conclusions. Two 
appendices deal with fundamental aspects concerning the conditions for PNC effects 
in the Floquet-eigenvalues. 



2 The Floquet operator 

In this section we will derive the Floquet operator which describes the time-evolution 
of the states of the subspace 71 for the case of external fields S, B periodic in time, for 
one oscillation period. The fields are furthermore assumed to be constant during finite 
time intervals and to change suddenly between them. 
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Let us assume a periodicity of the fields which sets in at a certain time to > 0: 

H,ff{to + nT + i) = H,ff{to + i), (2.1) 

with nGlN,0<t<T. The time evolution of some initial state |0) at time t = into 
\t) at time t > is governed by the Schrodinger equation (|1.3|) , whose solution is given 
as a time ordered exponential integral: 

\t) = r exp |-z dt'H^ffit')^ |0). (2.2) 

With respect to the basis (|1 . 1|) this may be written as a matrix equation, 

ia\t) = mt,0)U im (2.3) 

{a, (3 G {A, 5}; here and in the following we use the summation convention for repeated 
indices), with 

UM = ((7|rexp dt"H,ffit")^ |y)) , it > t'). (2.4) 

Due to the periodicity of H^^ff (pTTP, the evolution matrix U is also periodic: 

U{t,t') = U{t-T,t' -T), (2.5) 

hi to + T <t' <t <to + 2T, and thus 

U(to + nT + i, 0) = t/(to + i, to) Mto + T, to)]" Uito, 0) (2.6) 

(n G IN, < t < T). f/(to + r,to) is called the Floquet |T2|, matrix of the 
periodic Hamilton operator. It becomes immediately clear, that it is only necessary 
to know the time evolution of the state |0) for < t < to + T and the eigenvalues 
and corresponding eigenvectors of the Floquet matrix in order to calculate the time 
evolution up to arbitrary times. Indeed we have 

{uito + T,tor) =Y.i»\r)Kim, (2.7) 

r=± 

where A± are the eigenvalues of the (2x2)-matrix f/(to + T, to): 

A± = A ± v^, (2.8) 

with 

A = iTrI/(to + T,to), (2.9) 

K = A2-detf/(to + T,to). (2.10) 

Furthermore, {r\f3), (a|r), r G {+, — }, are the components of the left and right eigen- 
vectors of the in general non-hermitean matrix U(tQ + T, to): 

mto + T,to))^f,{f3\r) = Xr{a\r), (2.11) 
{r\a){U{to + T,to))a(, = ^r- (f|/3), (2.12) 



with the dual normahzation 

ir\s) = 6rs (2.13) 
and the completeness relation in the subspace TZ 

E = (2-14) 

r=± 

and with 

f/(to + T,to) = E IOAr(r|. (2.15) 

r=± 

Let us recall here that the diagonalization of a non-hermitean matrix U is always 
possible if its eigenvalues are all different [1^. If there is a degeneracy of eigenvalues 
one can in general only transform U to Jordan's normal form. In our case of (2x2)- 
matrices this reads 



(2.16) 



In the following we will in fact concentrate on the case where the Floquet matrix 
U{tQ + T, to) is of the form ( ^^.16] ), i.e. non-diagonalizable without the P-violating terms, 
but becomes diagonalizable when P-violation is included. This is precisely the case 
which is needed to get \/H^ P-violating effects. In summary: For the cases which we 
consider we always have non-degenerate eigenvalues of our Floquet matrix and thus 
the diagonalization (p.ll|) - (p.l5|) is guaranteed. For the points of special interest for 
us this non-degeneracy is only achieved when P-violation is taken into account. 

Let us now specify the configuration of electric and magnetic fields, which we will 
consider. We assume the fields to point in 3-direction and to be constant within certain 
time intervals. For this we choose times t^"^' (cr = 1, amax) with 

= < t'-'^ < ... < t*"'""^' = T. (2.17) 

Then, for < t < T, we write the electric and magnetic fields as 



max 



S{to + i) = ea ^ f(-'e(t'^^-t)e(t-t<'^-^0, (2-18) 

(T=l 

B{to + i) = ea^E S'-'e(t('^' -t)0(t-t''^-^O- (2-19) 



(7=1 



For arbitrary times t > to + T these fields are chosen to be 



£{t) = £{t-T), (2.20) 
B{t) = B{t-T). (2.21) 



Before to they are assumed to be constant: 



Sit) = ^'°'e3, (2.22) 
B{t) = iS'^'eg {0<t<to). (2.23) 



In Fig. we show an example configuration of {i^'""', 
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a=0 I 1 ' 2 ' 3 I 1 ' 2 ' 3 I 1 ' ^ 



to 







tp + r 

= T 



to + 2T 



Figure 1: Example of the time dependence of the electric and magnetic field strengths. An initial 
time interval < t < tg where the fields are constant is followed by an arbitrary number of Floquet 
cycles of period T. Each Floquet cycle is divided into amax (here: amax = 3) intervals of length 



To compute the time evolution of the system, we start with the effective Hamilton 
operator (|1.2|). In the subspace 71, w.r.t. the basis (|1.1|), it is represented by the 
(2x2)-matrix 

^max 

Heff{io + i) = E (2-24) 

where ( |2.24D defines H^.^ in the range to ^ ^ = + ^ < + ^- For later times it is 
defined by the periodicity relation ( |2.1| ). For < t < to, it is defined analogously, with 
the fields ( ^1221) , {^2^. 

In order to calculate the eigenvalues of the corresponding Floquet matrix we will 
perform the integration in the exponent of (2^) for each time interval with constant 
fields, where the time ordering plays no role, and then multiply the resulting matrices 
in the correct order: 

U{to + T,to) = r''l[ (2.26) 

(7=1 

where 

= exp {-tH^:J{&^ - &-'^)} . (2.27) 

THo- denotes the multiplication from right to left, with ascending a. Let us write the 
segmental (2x2)-matrices H^^^ (p.25|) using the SU(2) basis 

{(To= 12,^^;^ = 1,2,3}, (2.28) 
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where the cxj are the Pauh spin matrices, as follows: 



j=0 



■y) 



with 



Then (p.27|) may be expressed as 

3 

U^^^ = exp {-ihl^^At^^^} '^t^^i 

i=0 



with 



u. 



At'"-) 



cos 

-ih'f^ sin (ri'-^^A^^) /ri^^\ i = 1,2, 3, 

/ 3 \ V2 



(2.29) 



(2.30) 

(2.31) 
(2.32) 
(2.33) 



(2.34) 

(2.35) 
(2.36) 

(2.37) 
(2.38) 



Note that the 77'''' are even functions of H„. 



3 P-violating contributions to the Floquet eigen- 
values 

We would now like to ask the question, under what conditions there may be P-violating 



contributions to the Floquet eigenvalues X± (|2.8|). 

To study this consider the field configuration (|2.18| )- (|2.23|) and the one obtained 
by a reflection R on the 1-2-plane: 



•(a) 



(3.1) 



One signature of P- violation is 



AA± = A±({£:"'',B"''};if„)-A±({£:r,6r};i^w) 
^ 0. 



It is easy to see from ( p. 25]) that we have 



(3.2) 



(3.3) 
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and thus 

AA± = \±i{£'^\ B'^'}; H^) - \±{{£'-\ B'^'}; -H„). (3.4) 



Let us now consider the possibihty of AA± ^ 0. From (|2:26| ), (^^), (|2]2|) and 
(PI) we find 



det U{to + T, to) = exp | -2z ^ At*"' | 



(3.5) 



which is thus independent of if„ (cf. (|2.3CI|)). Therefore, taking into account 



( p.lOp , we only have to study the trace of U, as far as P- violating contributions are 
concerned. We will first show that for certain combinations of the electric and mag- 
netic field strengths' 3-components {£^''\ B'^'^^; a = 1, ...,ajnax}, there are no if„-linear 
contributions to the trace of the Floquet operator. For this we note that we can insert 
(2x2)-matrices V and V^^, with V non-singular, to the left and right of each t/'"', 
respectively. This leaves the trace invariant: 

r n W-^'] 



= Tr YI/"*"™""''!:^" • ■ • YfZ'''Y" j • (3.6) 

The matrices V and then act as a transformation rotating the Hamiltonians ^g^' 
in the exponents oiU^"^ (cf. ( ^.27| )). Consider now the case that the endpoints of all 2- 
dimensional vectors formed by the field strength components £''''\ B^"^ (a = 1, amax) 
lie on a straight line in the £^-;B-plane, i.e. they can be parametrized as follows: 



C\ (3.7) 



We show in Appendix ^ that in this case we can find a matrix V which, when applied 
as in ( p.6|) , transforms away all if^-odd terms in the trace of the Floquet operator. 

Thus we find that the "points" {£^^\B^'^^) must not lie on a straight line in the 
£^-i3-plane in order to obtain AA± ^ (cf. (p.4|)). This is, of course, only possible for 
a minimum of three time intervals, i.e. amax > 3. 

We will now provide an expression for Tr f/ ( pl6|) , for (Tmax = 3, in terms of h^^^ and 
the vectors u*"^', a = 1, 2, 3. Using the relation 

(aoO"o + a ■ cr){bo(TQ + b ■ cr) = (ao&o + a ■ b) ctq + (oob + boa + i[a x b]) ■ cr (3.8) 

we find: 

U{to + T, to) = W'WW =: E n,(l)a„ (3.9) 

i=0 

Uo{l) = |Tr (I/'^'f/'^'f/'^') 
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u(l) = i f/<^)f/«,cr} 



(!)/„, (3)^, (2) 



"0 



with 



(3.10) 

(3.11) 
(3.12) 



The arguments "1" refer to the number of Floquet cycles, which will be greater than 
1 later, where we will use again the representation ( |3.9| ) of the time evolution matrix 
in terms of the Pauli matrices. 

Taking into account (|2.30|) -( p33D and ( p.35| )- (|2.38| ) we find for the P-even and 
P-odd contributions Aq and Ai to A ( p.9|) the following expressions: 

(3.13) 
(3.14) 
(3.15) 



A = Ao + Ai 

Ao = I [A(/7„) + A(-/7„)] =: 

Ai = i [A(i7„) - A(-i7„)] =: Aie, 



with 



h,.,.h.,^to('>"'A'")„,(„..A,.)!"'('>"A"") 



(3) 



Ai = (h^^^ X h'^') ■ h 



+ cycl. perm, in 1,2,3 

sin(r7(^)At(^)) sin(r/(^' At'^') sin(r/(^)At(^)) 



(3) 



V 



(3) 



V 



(2) 



V 



(1) 



(h*^' X h'^') ■ h 



(3) 



■d{^A- /is) 



111 

gw £;{2) gi'i) 

^{1) ^(2) g{3) 



(3.16) 
(3.17) 

(3.18) 



For the eigenvalues A± of the Floquet matrix we now get from (p.8|)-( prT0| ), (p.5|), 
( CT ), and (|37[5| ): 



AA± 



e (^Ao + Ai ± VAg - 1 + 2AoAi + A? 
A ± v^, 

e (A^ - 1 + 2AoAi + A? 
e(2Ai± 



A^-l + 2AoAi + Ar'^' 



X2 
-^0 



l-2AoAi + At 



2^1/2- 



(3.19) 
(3.20) 



(3.21) 



We see that in general AA± will be linear in Ai, i.e. linear in 
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The basic idea which inspired us to investigate the questions presented in this paper 
was to explore the possibihty of sphttings AA± proportional to the square root of 
which, according to (|3.21|) requires 



X2 

-^0 



Ai ^ 0. 



(3.22) 
(3.23) 



This should yield P-violating effects larger by several orders of magnitude compared 
to ordinary if„-linear ones. As we will show this is indeed possible for unstable states 

A,B. It is not possible for the case of stable states (Ta = Tb = 0), for which the 

(<^) 



operators H^Zf are hermitean and U(to + T, tg) is unitary (cf. Appendix 



Next we will write down an expression for the iterated Floquet matrix U{to + T, to)^ 
(cf. ( |2.6|) ) in terms of the eigenvalues A±. Using the left and right eigenvectors (|2.11|) f. 
and their orthonormality ( |2.13D f. we find: 



U{to + T,toy 



For n = 1 this relation reads (cf. ( p.8|) ) 

Uito + T,to) = \l 



E \r)K{r\ 

r=± 

UK + ^-)'^ + UK - 



-)(- 



On the other hand we can express U in the basis 
with ( |3.25|) we find 



-)(+|-|-)(-| 



Hi- 

as shown in 



A 

l+)(+|-|-)(-| 

From (|3l9| ) we see that we can write A+ as 

A± = ^ exp(±£:), 



iTrI/(to + T,to) =Mo(l), 
1 3 

= Y^ Ui{l)(Ti. 
^ i=l 



tanhe: 



1 



Ao + Ai 



1/2 



A~' 



where e is in general complex. This gives 

A 



Thus we have 



= exp {±ne) = exp {—iuj±nT) 
u± = i{\n^±6) T-\ 



i=0 



(3.24) 

(3.25) 
Comparing 

(3.26) 
(3.27) 



(3.28) 



(3.29) 
(3.30) 
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with 

Uo{n) = Ccoshine), (3.31) 
u,{n) = Ui{l)Csinh{ne) /^/^, i = 1,2,3. (3.32) 

From (|3.3CI|) - (|3.32|) we see that for nearly degenerate eigenvalues (|k| <^ 1), there 
are oscillations of the matrix elements of the Floquet matrix where the (complex) 
frequencies uj± contain terms proportional to a/k. For fixed n and n 0, however, 
there are clearly no terms oc ^/K, in Ui{n) ( p.32| ), as we see from ( p.28| ). This is crucial 
for consistency: Consider the case ( p.22|) f., which we intend to study in the following 
sections. Then the oscillation frequencies become proportional to ^/K, ~ \/2Xi oc \/H^. 
But for fixed n, no odd powers of ^/H^ can appear in an expansion of U w.r.t. \/H^ 
since the Hamiltonian operator H^^ contains only a term linear in H„. 

We emphasize that \/H^ effects in the Floquet frequencies are only possible for 
unstable states. In Appendix we show that for stable states, i.e. for = = 
frequencies proportional to \/H^ are excluded. 



4 P-violating observables 

In this section we will study possibilities to obtain optimized observables to measure 
the P-violating contribution iJ„ to the Hamiltonian operator for systems of the form 
described in section ^. The general idea is to measure occupation numbers of the states 
\A) and \B) of the subspace TZ after a given number n of oscillations of the external 
fields. In other words, we will consider the moduli squared of some linear combination 
of the matrix elements of the evolution matrix U{tQ + nT,to). The initial state |0) is 
chosen to be either 1^4) or \B), and the fields S^^^e^, S'^'es, switched on between t = 
and t = to are used to provide a mixing of these states before the atoms enter the 
oscillatory field. In the following we will label this experimental setup as roman no. I. 
At t = to + nT the population of either 1^4) or \B) shall be measured. A measurement of 
such a population number for a sequence of n-values will reveal the Floquet frequencies 

Suppose we start with a state |/3) {j3 = A ot B) at time t = 0. We let this state 
evolve until a time to + , when the state vector of the undecayed atom is 

I/?, n- 1) = U{t, + uT, 0; B'^'}) |/5). (4.1) 

Its norm and thus the probability to find the atom undecayed is 

wi{f3,n) = {(3,n-\\f3,n-l). (4.2) 

The probability to find the atom in the state a {a E {A, B}) is 

wi{a,P,n) = \{a\l3,n-l)\\ (4.3) 

We also consider the evolution of the state \j3) in the reflected setup with the R- 



transformed fields £^^\B'--^^ (|3.1| ). The corresponding states and other quantities will 
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be labelled with a roman number 11 as opposed to those corresponding to the original 
setup, labelled by I: |/?,n;E), wn{P,n), etc. 

Suppose now that we repeat the experiment with the S'"^'} setup A^o,i times 

starting always at time t = with the atom in state \f3). The second experiment with 
the R-reflected setup {£^^\ ^r^} we repeat Nqj times. Let Nij{P,n) and Nij{a, P,n) 
be the number of undecayed and the number of atoms in state a (a G {A, B}), respec- 
tively, at time to + 

Let us assume that we choose the numbers of repetitions Nqj, Nqj such that the 
numbers of undecayed atoms at time to + riT are the same for each setup: 

Ni{l3,n) = N^{(3,n) = N{l3,n). (4.4) 

The relative occupation of |a) compared to the number of undecayed states, 
Nij{a, I3,n) /Nij{[3,n), is our observable. Its mean value is 

ri,n a, n) = — ,^ . — = ' ^ ^ ■ 4.5 

The expected numbers of repetitions necessary to achieve N{P, n) undecayed atoms 
are 

^N„,,} = ^^!^ (4.6) 

For an ensemble of N{(3,n) undecayed atoms the number of atoms in state \q) has a 
variance 

, 2 



(AAri,n(«,/?,n)) = (a^(/?, n) Ari,i(a, /?, n) 

= iV(/5,n)ri,n(a,/?,ri)(l -ri,i(a,/3,n)). (4.7) 

Suitable P-odd observables are the asymmetries between the relative occupation num- 
bers for the original and the reflected setups: 

v4(a, ^, n) = — — — — — — — . (4.8) 

ri{a,p,n) + rji{a,p,n) 



Let us now consider the statistical significance of the P- violating observable (|4.8| ). The 
variance of ( [4.8| ) for an ensemble of N{(3, n) undecayed atoms at time to + nT is given 
by 

4 rirn(ri + rn - 2rirn) 



(n,! = ro(a,/?,n)). (4.9) 



We can claim to see a l-a-effect if 



|A(a,/3,n)|^ > (Ay4(a,/5,n))^ (4.10) 
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Figure 2: Sketch of the real parts of the vectors h'"' for the choice given in Table I. 



Using ( [4 .61) this can be written as a condition on the mean total number of initial 
atoms in the two experiments with opposite chiralities of the field configuration which 
are necessary to establish a l-cr-effect: 



{No) = {No,i) + {Noj)> No 
No = 



1 ^ 1 \ 4riri(ri + rn - 2rirn) 



(4.11) 
(4.12) 



The task will now be to minimize the value of No as a function of n and the field 
strengths {S<-^\ B^^^; a = 0, a^ax}- 



5 Results 



In this section we present our numerical results. In the first place we calculate a solution 
for the fields {£^''\ B'"^ ; a = 1,2,3} for which the eigenvalues ( p.l9|) of the Floquet 
matrix satisfy the conditions ( p.22| ) f., i.e. are only split by a P- violating quantity 
proportional to \^H„. This corresponds to Ag = 1. Using this solution we calculate the 
asymmetry (|4.1|) and minimize the number A^o ( 4.12 ) of initial atoms required for a 1- 
a-effect, by choosing a convenient initial mixing through £^°\ B^°\ For our calculations 
we used the fixed constants (cf. |]1|) as given in (|1.7|). 



5.1 \/^-linear eigenvalues 



As can be seen from ( |3.16| ), ( p. IT] ), the conditions ( |3.22| ) f. may only be fulfilled for a 
convenient relative adjustment of the three vectors h*'^', a = 1,2,3 defined in (|2.29|) - 

(S). 
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(J 


1 


2 


3 


do / [KxlZ 1 




1 ^7 n 
iO ( .u 


S /17 1 n~2 
— 0.4/ • iU 


E^"'' [V/cm] 


-2.23 • 10^2 


41.3 


-2.23 • 10-2 


i[A + (^A-/UB)i3('^)]//i [kHz] 


157.1 





-157.1 


B^^^ [G] 


1.33 


1.48 


1.63 




1.59 


9.91 • 10-2 


1.59 



Table 1: Configuration of the electromagnetic field within the sections a — 1,2,3 of one Floquet 
cycle (cf. @f., (|23|)). At least one of the parameters has to be fine-tuned to more decimals than 
indicated in order to match exactly the degeneracy condition ( 3.20| ). The intervals Ai'"', a = 1,3, 
correspond to 2ttM'''> = 0.01 (kHz)"^ 



However, we will start by choosing appropriate time intervals At'*^' ( p.38|) . As it is 



clear from ( p.30| ) ff., ( [4.1| ) both the probabilities w±{l3^ n) ( [4.2| ) and the minimum initial 
atom number A^^o ( 4.12|) are governed by the overall factor ^, viz. are proportional to 



I^P'^ and 1^1 respectively. From ( |2T7| ), QOg ), (|3ll| ) we find 



lep" = exp{-i(r^ + rB)rn}, (5.1) 



3 



T = (5-2) 



(T=l 



In order to preserve the occupation numbers from being essentially zero before one 
period of the time evolution within the oscillatory fields, one has to choose the At*''' 
small enough. Stated differently, we would like to study the case, where a P-violating 
effect is accumulated during several up to many repeated cycles of the periodic fields. 
In this case the occupation numbers should show a P-violating oscillation with a wave 
number proportional to the imaginary part of As one derives from the numerical 
values given in (p..7[), a Floquet time T (|5.2| ) of 



T = 3.5/is (5.3) 

yields a decay by a factor of |.^|2" ^ 1/10 within n = 10 periods. 

In order to have a non-vanishing Ai (|3.17|) one has to choose the fields and time 



intervals At'''' in a way, that the sines in (|3.17| ) are non-zero. On the other hand we 



would like Aq to equal plus or minus 1. As one finds numerically, for the total sum of 
the time intervals At'"', a = 1,2, 3, not to exceed the value ( ^.3|) , this is only possible, 
if the sines as well as the cosine in the second term of ( p.l6| ) are close to +1 or —1. 



(Then the two other terms of the sum over the cyclic permutations in 1, 2, 3 in (|3.16|) 
are close to zero.) A different choice, where more than one cosine cos(?7''''At'''') is away 
from zero, leads to a reduced, non-optimal value of the modulus of Ai. Now, if two of 



the sines are to be close to 1 and (|573|) is supposed to be the upper limit of (|5.2|) , the 
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Figure 3: The probabilities wi,i(a, f3, n) (4.3) ior a — B, (5 = A, H^/h — 1 kHz, as a function of the 
number of Floquet cycles n. The electric and magnetic fields are chosen as in Table |. No premixmg is 
chosen: to — 0. The open boxes indicate the values of wi(a, /?, n), the filled boxes wji{a, 13, n), i.e. those 
for the R-reflected setup. The dashed hues represent the complete time evolution of |(a|C/(t, 0)|/3)p 
as a function of t which is given in units of T. This plot illustrates the first Floquet beat for an 
hypothetical value for H„, where the overall decay does not yet overwhelm the P-violating effect in 
the first beat. 



electric and magnetic field terms dS'^'"\ [A + (/i^ — /is)i3''^']/2 in the corresponding two 
segments have to be chosen large compared to the imaginary part (Ta — Tb)/4: of h'^^'' 

id. (pi]), (Pi, (IOTP ). 



We have chosen a configuration as given in Table [^. The representation of the 
projection of the real parts of the vectors h'"^' is moreover sketched in Fig. p. From this 
it becomes clear, that the dot product h ■ h of the "unit" vectors h := h*'^'///''^'' 
in the relevant term of ( |3.16D is close to —1 (only close, because the imaginary parts of 

^ (1) (2) ^ (2) ^ (3) 

the introduce small perturbations). Similarly we have h ■ h ~ 0, h ■ h ~ 0. 
Then we get for Aq from ( |3.16| ): 

Ao ^ cos{r]^'^At^'^) cos{r]^'^At^'^) cos{r]^'^At^'^) 

+ sm{r]^'^At^'^)cos{r]^'^At^'^)sm{r]^'^At^'^). (5.4) 

Now we make the factor cos(r7^^^At*^^) to be close to 1 by choosing the time interval 
At*^' small compared to the two other ones. Moreover we choose r^^^'At'^' ^ r^'^^'At'^' 
(cf. Table p. In this way we clearly have Aq — 1 and by fine tuning of the parameters 
we can achieve Aq = 1- This choice of parameters yields that sin(?7^^'At^^')/r7^^' ~ At'^* 
and hence that Ai is proportional to At^^K We have not found a way to avoid such a 
relatively small proportionality factor, since the alternative choice of a small 77'^' would 
result in an equivalent reduction in the triple product of the vectors h*"^^ in Ai, which 
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Figure 4: The logarithm of the probabiUties 'Wij{a, f3, n) (4.3) for the same configuration as in Fig. 
^, as a function of the number of Floquet cycles n. Compared to Fig. the first four Floquet cycles 
are shown. The oscillatory, i.e. the lower sequence of points refers to wji{a, P,n), the upper one to 

wi(a, /?, n). 



is equal to the area of the triangle spanned by the projections of these vectors on to 
the 1-3-plane (cf. Fig. ^). 

In summary, using the values given in Tab. we obtain: 



Ao 
Ai 

UTa + Tb)T 



1.0, 
1.25 
2.49 
4.46 



10" 
10" 
10" 



(5.25^0.108) kHz, 
(5.25^*0.108) kHz, 
2.17- 10"^ 



(5.5) 

(5.6) 
(5.7) 
(5.8) 
(5.9) 
(5.10) 



(Here Ai is real and positive by chance, which is not true for a general choice of the field 
configuration.) As we have mentioned before the fraction of the P- violating oscillation 
wave number ( ^.7| ) divided by the total decay width (|5.10|) is 



70. 



4e 



{Ta + Tb)T 



2.06-10" 



(5.11) 



Our numerical investigations showed that this cannot be increased by more than a few 
tenths of a percent. It is essentially governed by the fraction ^/H^/(TA + F^). Longer 
lived states with smaller decay widths will not a priori lead to a larger fraction ( ^.11|) , 



since then At*^^ has to be smaller in order to fulfill ( |5.5|) . 



16 



0.6 



0.4 - 




d -0.2 - 



-0.4 - 
-0.6 - 
-0.8 - 

-5 5 10 

d^cy/i [kHz] 

Figure 5: The asymmetry A{a,P,n) for a = /3 = A, 8 < n < 12, {5'"' , B'"', Ai'"'; cr = 1,2,3} 

as given m Tab.|l|, A + (/xa-mb)6<'" =0, as a function of d£''°^ . Contributions of order have been 
neglected. The curves are essentially zero except for the small intervals, where resonances occur. The 
approximate centres of these intervals are given in Table |[ Each resonance represents a different n. 
Hence, the largest one at dS'-"^ /h ~ 1 kHz corresponds to n = 10. 



For illustration we show in Figs. ^, H the probabilities w\^]^{a, /5, n) ( [4 .31) for a = B, 
j3 = A as a. function of the number of Floquet cycles n for an hypothetical P- violating 
matrix element of H^/h = IkHz. In this case, the overall decay determined by ^ is 
weak enough such that the first Floquet beat shows a clear P-violating change under 
a reversal of the fields' handedness. In Fig. |^ also the complete time evolution of the 
probabilities |(a|Lf(t, 0)|/3)p (cf. ( p.4|) ) in between the Floquet points is shown. Fig. | 
illustrates that the oscillations may occur for only one chirality (cf. ( |5.8| ), ( p.9|) ). 

In summary, with the arrangement investigated so far it will not be possible to 
actually observe P-violating oscillations of the occupation numbers, which here, due to 
the reality of Ai, would only show up for one handedness of the fields. On the contrary, 
most of the atoms will have decayed long before one such oscillation is completed. 
However, as we will show in the following, it should be possible nevertheless to observe 
large P-violating effects. 

5.2 Results for A{a,/3,n) and Nq 

In this section we give the results of an optimization procedure for the asymmetry 
A{a,(3,n) ( [4.8|) and the minimum initial number Nq ( [4.12|) of atoms required for a 
l-(T-effect. For this we introduce a "premixing" of the states between t = and t = to 
by applying electric and magnetic fields S^'^^e^, B^^^e^ (|2.22|) . 

We considered an evolution of n < 20 Floquet cycles, which corresponds to a time of 
approximately 60 /xs (cf. ( ^.21) , Table |l]). Within this period, 98% of the initial number 
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Figure 6: The decadic logarithm of the minimum number Nq (4.12) of initial atoms for \a) = |/?) = 
\A), n = 2,6, 10, 14, 18, for the same configuration of electric and magnetic fields as in Fig. ||. The 
scale of the horiz ontal axis is stretched and denotes the difference dAf'"' from the shifted central 
values df'"' (Isl^ ) as given in Table || for each n. Contributions of order have been neglected. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


dSlP^ /h [kHz] 


67.7 


61.2 


53.7 


45.1 


36.1 


27.4 


19.0 


11.9 


6.1 


1.1 


n 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


dSjfyh [kHz] 


-3.0 


-6.4 


-9.1 


-11.5 


-13.5 


-15.2 


-16.7 


-18.0 


-19.1 


-20.0 



Table 2: The central values d£^°'' (cf. ( 5.12 )) of the resonances in Fig. I for n= 1,...,20. 



of states have decayed. However, for n < 20, the modulus of the P-violating arguments 
of the hyperbolic sine and cosine in ( |3.31|) f. is smaller than 0.05. In Fig. ^ we show the 
asymmetry A{a, f3, n) (g^) for a = /? = A, 1 < n < 20, {£^^\ B'^^\ At^''^; a = 1,2, 3} as 
given in Table [|, choosing S"" = -A/{ha-I^b) = 1.4762 G, i.e. A + (/iA-/iij)i3<'" = 0, 
as a function of dS'^°\ In computing A{a,l3,n) we neglected for the first contributions 
of order H^. We find an optimal asymmetry for |a) = = \A), n = 10, and 
~ 0.28 V/cm of A{a, (3, n) = -0.84. 

In Fig. 1^ we show the decadic logarithm of the minimum number Nq ( |4.12| ) of initial 
atoms for \a) = = \A), n = 2,6, 10, 14, 18, for the same configuration of electric and 
magnetic fields as in Fig. |^. Note however, that for an appropriate resolution the scale 
of the horizontal axis is stretched and denotes the difference dAS^°^ from the central 
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Figure 7: The decadic logarithm of the fraction A^o/|^(a, /?, n)| of initial atoms for \a) = \(3) = \A), 
n = 2,6, 10, 14, 18, for the same configuration of electric and magnetic fields as in Fig. ||. The scale of 
t he h orizontal axis is stretched and denotes the difference dAE'-"'' from the shifted central values d£^"^ 
(5.12) as given in Table for each n. Contributions of order have been neglected. 



values d£^J'\ shifted by (10 - n)/4kHz/i: 

dS'-"^ = d£f^ + ^(10 - ^)kHz h + dA8^'\ (5.12) 

as given in Tab. ^ for each n. 

In Fig. 1^ the decadic logarithm of the fraction NQ/\A{a, I3,n)\ is plotted for the 
same parameters as in Fig. ^, again as a function of A£^'"'. We find that, although the 
number of atoms Nq increases when n becomes larger, as can be expected from ([4.12|) 
which is proportional to the factor (cf. ( ^.11) , ( ^.101) ) 

ler'" = exp{0.22 n} = lO°-094-^ (5.I3) 

the observable A{a,l3,n) at n = 10 is still the optimal one in the following sense: Al- 
though the asymmetry at larger values of n may have a higher statistical significance 
indicated by a smaller value of Nq, the smaller value of the observable might be dom- 
inated by systematic experimental errors. In this sense, the fraction of A'o divided by 
A{a,l3,n) may provide a useful measure for an optimal compromise when taking into 
account different sources of errors. 

In Figs. |a-d we plot A{a, (3, n) (gj), Ig JVq, lg(]Vo«;i,n(a, /3, n)), \g{No/\A{a, /3, n)|) 
for n = 10, a = (3 = A, and the other parameters chosen as in Figs. ^|-^ From Fig. §3 
we can see that dS^^^h = 1.07kHz (f^*"^ = 0.28 V/cm) gives the minimal value for the 
initial number of atoms required for a 1-cr-effect: 

]Vo = 4-10^ (5.14) 
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Figure 8: (a) A{a,P,n), (b) IgA^o, (c) \g{Nowij{a, P,n)) and (d) lg{No/\A{a, (3,n)\) for n = 10, 
a = P = A, and the other parameters chosen as in Figs. ||-|^. In (a) the dashed hne represents 
the approximation, where terms of order are neglected (cf. Fig. |^), the solid hne is the exact 
asymmetry (18). In Fig. (c) the solid hne represents lg(7VoU'i), the dashed hne \g{NoW]i). Thus (c) 
shows how the occupation number changes under a reversal of the fields' handedness. 



There we have an asymmetry of 

A{a, P, 10) = -0.7. (5.15) 

Fig. ^ shows that, starting with 4 ■ 10® atoms, after n = 10 steps, one should observe 
either about 10 or only 2 atoms, depending on the handedness of the field configuration. 
Fig. ||d shows that lg{NQ/\A{a, P,n)\) is also almost at its minimum. 
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Figure 9: The asymmetry A(a,/3,n) (gj) for a = /3 = A, n = 10, B*"', At'"'; cr = 1,2,3} as 

ffiven m Tab. g, A + {^lA - ^iB)B^°^ = 0, as a function of d£'°^ for different values of the P-violatinp; 
matrix element H„. 



In Fig. ^ we plot the asymmetry A{a,P,n) (|4.8| ) for a = P = A, n = 10, and the 
other parameters as in Fig. |^a, for different values of the P- violating matrix element H„ . 
This illustrates, that H„ can be determined from the magnitude of the asymmetry and 
the positions of its extrema. We point out that also the sign of H„ may be determined 
from the asymmetry, since it changes sign under H„ — ^ — -ff„. 

We also varied the magnetic field and found that slightly away from the Zeeman- 
level degeneracy A + (//^ — = one can even reach a 100%-effect! However, the 
number of atoms required to observe a 1-cr-effect is of the same order of magnitude as 
above. Thus we will not investigate this further here. But it shows that the occurence 
of large P- violating asymmetries is not necessarily connected to the P-even degeneracy 
of the Floquet eigenvalues. 

To summarize, our results of this section show, that for a choice of parameters as 
given in Tab. ^ large asymmetries may be reached easily within a few Floquet cycles, 
for premixing fields within an accessible range. 
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6 Conclusions 



In this article we have investigated parity violating effects in atomic dysprosium where 
one has a near degeneracy of two levels of opposite parity. We considered the time evo- 
lution of this 2-level-system in external electric and magnetic fields having a periodic 
structure in time with period T. The object to study was the Floquct matrix which 
describes this periodic time evolution. We found that for unstable states the eigen- 
values of the Floquet matrix which give the eigenfrequencies of the system can have 
contributions proportional to the square root of the P-violating weak interaction matrix 
element H^. This leads to beat frequencies proportional to \/H„. For dysprosium in 
the simple arrangements of external fields which we considered we found, however, that 
it is difficult to observe such beat effects since the states decay too fast. On the other 
hand, we found that very large P-violating asymmetries are easily obtained after a few 
Floquet cycles. This is achieved by tuning carefully the initial electric and magnetic 
fields in which the system evolves before entering the Floquet cycles. The electric and 
magnetic field, the time constants etc. necessary for this are all in an experimentally 
accessible range. The number of atoms necessary to obtain a statistically significant 
signal is generally of order 10^ to 10^ if \H„\/h ~ 2 Hz. Once having measured the 
asymmetry one may determine both the modulus and the sign of H„. It is clear that 
our considerations can easily be extended to other periodic electric and magnetic field 
arrangements and other 2-level atomic systems. We hope that these ideas will stimu- 
late experimental work in this direction. 

Acknowledgements: The authors are grateful to M. Zolotorev and D. Budker for 
discussions and correspondence which led to these investigations. They would also like 
to thank T. Brunne for useful discussions and suggestions. 
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Figure 10: Sketch of the vectors (f-"' , S'"') for different tr, chosen such that they he on a straight hne 
in the S-B-plane. For such configurations the Floquet eigenvalues have no P-violating contributions. 



Appendix 

A Explicit form of the matrix V 

In this appendix we give the exphcit form of the matrix V (cf. ( |3.6|) ), which rotates the 
matrices U'-''\ a = 1, (Jmaxi such that the trace ( |3.6|) looses any /7„-odd contributions 
if ( |3.7| ) is satisfied. 

We start with a set of electric and magnetic field strengths' components 
{£^*'^', S'"'; a = l,...,amax} lying on a straight line in the £^-i3-plane (cf. ( |3.7| ) and 
Fig. |10D. Of course, vs,b^'U)s,b and C'"^' in (|3?^ ) are real. To have a nontrivial field 
configuration we require 

wl + wl^Q. (A.l) 
The corresponding vectors h'"^^ (( |2.31| )-( |2.33| )) are then 

+ v + wC^^\ (A.2) 

dw£ \ 
. (A.3) 

-i(/iA - ^■b)wb ) 

Without loss of generality we may assume 

w ■ w = 1, V ■ w = 0. (A. 4) 
The sought matrix V may be composed out of two separate transformations, 

V = V2Y^. (A.5) 
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We choose the first rotation to be one about the 2-axis, such that w becomes parallel 
to the 1-axis. The second one is chosen to be a rotation about the 1-axis. This will 
bring all vectors ([A .21) into the 1-3-plane. Our ansatz thus reads: 



Choosing a such that 
we find: 



exp |fa(T2} = 
exp {f/3o-i} = 

cos a = Wi, 



cos(a/2) sin(a/2) 
— sin(a/2) cos(a/2) 

cos(/3/2) isin(/3/2) 
isin(/3/2) cos(/3/2) 



sm a 



-W3 



i=0 



with 



hi 



''■0 5 



h'l' cos a — h'-^^ sin a 



h'r 



-W3 

h':^^ COS a + sin a 



2 ' 4 



-Wi 



iA + Kr^-r^) +(vxw)-e2=:/i^. 



(A.6) 
(A.7) 

(A.8) 
(A.9) 

(A.IO) 

(A.ll) 
(A.12) 

(A.13) 



Note that h2''^ and h^"^ are the same for all a. The angle (3 of the second rotation is 
chosen to satisfy: 



cos (3 
sin 13 



(h'sY + Hi 

(Kr + Hi 



■1/2 

-1/2 



(A. 14) 



Here we assume h'^ 7^ -^iH„ which will hold except for very special cases. It is clear 
that (3 is in general complex, near zero and odd under the exchange H„ —>■ —H„. Using 
O, O) and (lAlgp- dAl^) we find 



with 



-Y2 






\i=0 


'If \ 


— '''0 




= h< 




= 0, 


Hi \ 

h,<" 





iKr + Hl 



(A.15) 



(A.16) 
(A.17) 
(A.18) 
(A. 19) 
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We see from ( [A.16| )-( |A.19D that the /i/"^', j = 0, ...,3, are all even functions of H„. 
Inserting hj^"^ instead of hj'^ in ( |2.29| )-( |2.38| ) we see that also the corresponding Floquet 
matrix U"{to + T, tg) is an even function of H„. The same is then true for its trace 



which is the same as the trace of the original Floquet matrix according to (|3.6| ) 



Tr [/(to + T,to) = TrU"{to + T, to) = even function of 



(A.20) 



Using also (|33| ) and (|2.8| ) - (|2.1CI|) we prove our assertions made in Section ^, that if the 
points {S^^\B^'^^) lie on a straight line, the Floquet eigenvalues are even functions of 
H„ and thus have no P-odd contributions. 



B The case of stable states, = = 

In this appendix we show explicitly for the case of a time-evolution with amax = 3 
sections of constant fields within one period of oscillation, that for stable states, i.e. 
Ta = Tb = 0, conditions ( p.22| ) and ( p.23| ) may not be fulfilled simultaneously. As 
stated in section ^, Ai ( p.l7|) will vanish in any case where Aq (p.l6|) is equal to 1. 

The reason for this becomes clear from the following general argument: For stable 
states the Floquet matrix U is unitary. A unitary matrix U may always be written 
as U = exp{—iHT} , where H is a hermitean matrix. For the case that U and H are 
(2x2)-matrices, the eigenvalues of H are given as 



uj± = lTrH±^J{lTrHy-detH. (B.l) 

Let the series expansion of matrix elements of H'-°^ w.r.t. the small parameter 6 := 
H„/A, where A is some constant with the dimensions of be Hij = H-j^ + SH^j'' + 
0{6'^). Then the expansion of the eigenvalues ( p.l|) reads 



u^± = \{Hfl + H^)±\{{Hfl-H^)' + AHf^H, 



(0) 

21 



+ 25 



1/2 



+ 0{5). (B.2) 

Now, if the eigenvalues are degenerate to the order 5°, we may diagonalize H such that 
Hfi = H22 and -^{2' = -^2? = 0- Then it becomes immediately clear from ( p.2|) , that 
the radicand is zero up to order 5. Hence, no contributions to uj± proportional to \/5 
are possible. The same is then true for the eigenvalues of U. Note also that, from the 
physical point of view, such contributions would make uj± complex for one handedness 
of the electric and magnetic fields, which can not be true since H is hermitean. 

We now want to prove that for = = the expressions ( p. 161) and (|3.15| ) have 



this property for any field configuration. In order to do this we will first verify a lemma 
for the trace of the product of two unitary matrices f/''^-' of the form ( p.34| ). Let us first 
introduce the notation 

s^"^ := sin(r/^"^At("^), 



25 



(^) 



cos(r/<"'At<"'), 

H — 2 — , 



cos 



h 



- (<t) - (r) 

h ■ n . 



(B.3) 



Note that for stable states, t]'-"^ ( p. 371 ) is the length of the real vector h''^' and thus h 
is a unit vector and c*''^' the cosine of the angle between two such vectors. We will 
prove now the 

Lemma: Given the half of the trace of the product of two unitary matrices U'-^\U'-^\ 
with determinants detU^^^ = l,a = 1,2, i.e. = (cf. (|2:29|)-(g38D), 



T, 



21 



iTr iW'^U 



its modulus may then and only then equal one. 
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if either 



.(2)1 



.(1)1 



(B.4) 

(B.5) 
(B.6) 



or 



.(21) I 



1, (B.7) 

i.e. if both matrices equal plus or minus the unit matrix or h*^' and h*^' are parallel to 
each other. 

For the proof one first calculates 



1 - c'^'c'^' 
1 + c'^'c'^' 



3(+)2 _ ^(-)2 



From this we obtain the equivalences: 



21 



21 



-1 



1 - d'>c' 



(2)^(1) 



(-)^ 



1 + c^'^d 



+ s 

(2)„(1) 



.(21) / „(+) 



' — s 



(-)' 



.(+)^ 



+ c 



(-)" 



_ ^(21) |^^(+)2 . 



.(-)' 



(B.8) 
(B.9) 
(B.IO) 
(B.ll) 

(B.12) 
(B.13) 
(B.14) 
(B.15) 



In the case that |c^^^^| < 1, it follows immediately from both ( |B.13|) and (|B.15|) — since 
< < 1 and < c'*'^ < 1 — , together with (pAO| ), (IRUD , that 







(B.16) 



26 



and thus from ( P?!^ and that (|B;B|) is necessary for {^). For c''^' = ±1 ( [KTI) 

it follows from (pTT^ ) and (|B7[5|) that Tis = 1 for s^^^ = and T21 = -1 for c'^' = 0. 
qed. 

By applying this Lemma to the trace of the product of a third matrix f/'^* with 
the product of the first two matrices, we will now prove the initial statement. In the 
follwing H„ will be set to zero if not otherwise stated. Then the (real) vectors h*"^', 
(7 = 1, ...,3, all lie within the 1-3-plane (cf. ( |2.31| )-( |2.33| )). Consider the product 



■ (2) 



■ (2) " (1) 



.(2)^(1) 



(7. 



Defining 

one finds (cf. ([B^)) 

|h(2l)|2 = I-T2V 

Therefore ( [B.17| ) may be written as 



f/(2i) = ]lcos(|h'^^)|) -2h*''^cr sin(|h(^i)|) 



(B.17) 

(B.18) 
(B.19) 

(B.20) 



■ (21) 



with h = h(^^V|h(^^)|. 

Let us now consider the trace (divided by 2) 



321 = |Tr (f/'^'f/(^'f/' 
= |Tr (f/'^'f/(^^)) 



(B.21) 



and investigate two cases: (a) = 1, i.e. s'^^' = (1 — c^^^'^)^/^ = 0, and (b) |c*^^'| < 1 
i.e. s(^i) ^ 0. 

In case (a) we have 

h(21) (C(2)^(l)^3(2)^(l))_ 

From the above Lemma we therefore conclude that 



is only possible if either 



321 1 



(3) (21) , 



h h 



and thus all "points" h*"^', cr = 1, 2, 3, lie on a straight line, or if 



.(3)1 



cos(|h 



(21)1 
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(B.22) 

(B.23) 
(B.24) 

(B.25) 



and thus s*^' = 0. Note that ( |B.25D may be fulfilled again due to the Lemma. 

In case (b) the vector h*^^' does not lie in the 1-3-plane which is spanned by the 
vectors h'^' and h'^' because of the cross product term in ( B.18|) . Then, again due to 
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the above Lemma, ( |B.23| ) requires ( |B.25[ ) to be fulfilled since ( |B.24D is not possible as 
h'^^ lies within the 1-3-plane (remember that = 0). Then s'^' = s'^' = due to the 
Lemma, and s^^' = from (p.25| ). 

In summary, in any case where ( [B.23| ) is fulfilled for stable states, i.e. Ag = 1 (cf. 
( P.16| ), (|3.19|) ff.), the if„-linear term 2^^AoAi, and thus e (|3.28| ) vanishes as well, either 
by the arguments presented in Appendix |^ or by the sine s*^' being zero (cf. ( |3.17| )). 
qed. 
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